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ABSTRACT 

One of the main astrophysical processes leading to strong emission of gravitational waves to be detected by 
the future space-borne interferometer LISA is the capture of a compact star by a black hole with a mass of a 
few million solar masses in the center of a galaxy. In previous studies, main sequence stars were thought not to 
contribute because they suffer from early tidal disruption. Here we show that, according to our simulations of 
the stellar dynamics of the Sgr A* cluster, there must be one to a few low-mass main sequence stars sufficiently 
bound to the central Galactic black hole to be conspicuous sources in LISA observations. The probability that a 
white dwarf may be detectable is lower than 0.5 and, in spite of mass segregation, detection of a captured neutron 
star or stellar black hole in the center of the Milky Way is highly unlikely. 

Subject headings: black hole physics - Galaxy: nucleus - gravitational waves - stellar dynamics 


1. INTRODUCTION 

Recently, the presence of massive black holes (MBHs, with 
masses from a few 10 ^ to a few IO^Mq) in the center of non¬ 
active galaxies has received considerable observation support, 
mainly from the kinematics of stars or gas in the center-most re¬ 
gion of nearby galaxies (Kormendy & Gebhardt 2001). One of 
the most convincing cases is the Milky Way for which the mo¬ 
tion of individual stars can be observed, both along the line of 
sight or perpendicular to it (Genzel et al. 2000; Ghez et al. 2000; 
Schodel et al. 2002). They orbit a dark mass of 2.6 x IO^Mq 
whose extent is less than 0.001 pc and coincides with the radio 
source Sgr A*. 

These MBHs reside at the center of stellar clusters whose 
densities may exceed IO’MqPC"^ (Lauer et al. 1998; Alexan¬ 
der 1999). An exciting consequence of this is the possibility for 
a star to be captured onto a relativistic orbit around the MBH by 
emission of gravitational radiation (Hils & Bender 1995; Sig- 
urdsson&Rees 1997; Freitag 2001; Ivanov 2002). Forasource 
at a distance of a few hundreds of Mpc and a MBH’s mass in the 
range between 10 ^ and IO^Mq, the frequency of these waves 
will enter LISA’s range as the orbit shrinks down (Thorne 1998; 
Danzmann 2000). If it withstands the tidal forces, the star will 
eventually find itself on an unstable, plunging orbit and disap¬ 
pear through the BH’s horizon. The computation of the late, 
strong-held, phases of this inspiral requires full general rela¬ 
tivistic treatment and cannot be handled with tools presently 
available (see Glampedakis & Kennehck 2002; Glampedakis 
et al. 2002 and references therein). Nevertheless when the 
pericenter distance, Rp is still well beyond the stability limit 
(3Rs = 6GMbhc~^ for a circular orbit around a non-rotating BH, 
with G the gravitational constant, Mbh the BH’s mass and c the 
speed of light), approximate methods can be used to determine 
how the orbit evolves through emission of gravitational waves. 
For a star with mass M* on an initial orbit of semi-major axis 
flo and high eccentricity, 1 -eo “C 1 , the time to plunge is ap¬ 
proximately (Peters 1964) 
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where is the initial pericenter distance. Compact stars, 
i.e. white dwarfs (WDs), neutron stars (NSs) or stellar black 
holes (SBHs) spiral all the way down to the horizon but a 
main sequence star (MSS) is torn apart by the tidal gravita¬ 
tional field of the MBH at a distance Rj = A:(Mbh/Mh<)'/^Rh< ~ 
5OkRs(MBH/lO6M0)-2/3(p,/l gem ^) where R* is the ra¬ 
dius of the star, p* its average density and k, a constant of order 
unity, depends on the stellar structure (Hills 1975; Rees 1988). 
Hence a MSS may come close enough to the MBH to emit sig¬ 
nificant amounts of gravitation radiation only if it is compact 
enough. Note that p* is maximum for M* ~ 0.07 Mq, at the 
transition to brown dwarfs (Chabrier & Baraffe 2000). 

We only consider quadrupolar gravitational radiation be¬ 
cause it dominates energy and angular momentum losses. In the 
weak field approximation and neglecting the gravitational influ¬ 
ence of other stars, the orbit may be treated as a Keplerian el¬ 
lipse whose parameters slowly change as a result of emission of 
gravitational radiation. The latter is emitted at integer multiples 
of the orbital frequency, u}„ = ^(GMbh)*^^^”^^^- At a distance 
D from the source, the strain amplitude in the n* harmonic is 
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where the non-dimensional factor 7 is an intricate function of 
n and the eccentricity e (Peters & Mathews 1963; Pierro et al. 

2001) . For simplicity, we consider the the nns amplitude, aver¬ 
aged over-l- and x polarizations and all directions. 

We assume that a star is irremediably captured by the MBH 
if it gets on an orbit with a time scale for shrinkage by emis¬ 
sion of gravitational radiation shorter than the time over which 
2 -body relaxation could significantly modify the pericenter dis¬ 
tance, Tmod — 0^7iei < Tiei (Sigurdsson & Rees 1997; Ivanov 

2002) . Trei is the usual relaxation time (Binney & Tremaine 
1987) and 9 is the angle between the trajectory and the direc¬ 
tion to the center. The typical eccentricity of capture orbits may 
be estimated from the condition Tpiunge < TreP 
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Here, Vorb is a typical velocity for stars in the “cusp”, i.e. the 
central region of the cluster where the kinematics is dominated 
by the gravitational attraction of the MBH. Mcusp, of order Mbh, 
is the total stellar mass in this region. The rate of capture is con¬ 
trolled by processes that send stars onto these extremely elon¬ 
gated orbits. In a spherical cluster, the most important one is 
2 -body relaxation. 

2. NUMERICAL MODEL OF THE CENTRAL CLUSTER 

To compute capture rates, many difficulties have to be faced, 
in particular the role of mass segregation, stellar collisions or 
tidal disruptions. Although these complications can be included 
to some extent into analytical computations (Sigurdsson & Rees 
1997; Miralda-Escude & Gould 2000), the intricate nature of 
the problem calls for numerical simulations. We have recently 
developed a new Monte Carlo (MC) code to follow the long 
term evolution of galactic nuclei (Freitag & Benz 2001, 2002; 
Freitag 2001). This tool is based on the scheme proposed by 
Henon (1973) to simulate globular clusters but, in addition to 
relaxation, it also includes collisions, tidal disruptions, stellar 
evolution and captures. 

For this work, we devised a model which represents the cen¬ 
tral cluster of our Galaxy. It is an ry-model (Tremaine et al. 
1994). The stellar density is p^,(R) = r]Mci/(4-7TRl)r'^~^(l + r)^~^, 
with r = R/Rb, r] = 1.3, Md = 8.67 x 10^Mq (total stellar 
mass) and Rb = 22 pc. The central BH has a mass of Mbh = 
2.6 X 10® Mq. These quantities were chosen to provide a good 
agreement with the observed run of the enclosed mass around 
Sgr A* as a function of the distance and the star counts (Genzel 

def 

et al. 2000), while imposing p. = Mbh/Mci = 0.03, a relatively 
large ratio chosen for the sake of resolution. We assume all 
stars formed 10 ^° years ago with a “universal” initial mass func¬ 
tion (IMF), chV*/dM* oc M““ with a = 0.3 between 0.01 and 
O.O 8 M 0 , 1.3 between 0.08 and Q.SMq and 2.3 up to 120 Mq 
(K roupa 2001). We do not include giant stars. All WDs, NSs 
and SBHs are assumed to have 0.6, 1.4 and 7 Mq, respectively. 
There is no initial mass segregation. The number of particles 
for our main simulation is 6 x 10 ®, so that each particle repre¬ 
sents 65.5 stars. 

3. RESULTS OE SIMULATIONS 

In Fig. la, we show the capture rates as a function of time. 
As captures are rare, we have integrated the evolution of the 
nucleus over many billion years to improve statistics. However, 
over such a long time scale, the nucleus model experiences a 
relatively important evolution. Most notably, significant mass 
segregation occurs. At the end of the simulation, the density 
is dominated by SBHs (which represent only 2 % of the stellar 
mass) in the region interior to ~ 0.15pc while the density of 
MSSs has decreased from a steep p* oc to a milder 
peak. This explains the steady decrease in the MSS capture 
rate. Given its efficiency, mass segregation should in princi¬ 
ple be introduced from the beginning of the simulation, which 
is unfortunately impossible for lack of observational constrains 
and adequate cluster models. Anyway, it is clear that the pre¬ 
dicted MSS capture rate is in excess of 2 x 10“® yr“^ and may be 
as high as a few 10“®yr“'. For WDs, the rate is a few 10“^yr“^ 
and around 5 x 10“* yr“^ for NSs and SBHs. These latter values 
are based on a small number of events and have low statistical 
significance. In Fig. 1(b), the initial orbital parameters are re¬ 
ported for all captures. One sees that orbits are very elongated. 


as predicted, with very small pericentre distances. Previous, an¬ 
alytical, studies did not address the capture of MSSs. They pre¬ 
dicted a WD capture rate ranging from 10“* yr“' (Hils & Bender 
1995; Ivanov 2002) to 10“^yr“^ (Sigurdsson & Rees 1997). Ex¬ 
plaining why these estimates are different from each other and 
from our result is beyond the scope of this letter as it would re¬ 
quire a lengthy discussion of the different nucleus models and 
various treatments of the physics in these papers. We note that 
our WD rate is similar but larger than that of Sigurdsson & Rees 
(1997) while these authors predict an initial SBH capture rate 
some three orders of magnitude larger than ours! 

During the course of the MC simulations, captured stars are 
immediately removed from the cluster. Their initial orbital pa¬ 
rameters being known, we come back to this data and compute 
the inspiralling and concomitant gravitational waves emission 
(Pierro et al. 2001; Glampedakis et al. 2002). To illustrate this, 
in Pig. 2, we plot the evolution of the frequency and amplitude 
of the lowest five harmonics of the emitted waves for typical 
WD and MSS events. The WD would be detectable during a 
few 10® years and the MSS for nearly 10® years. NSs and a 
SBHs typically have shorter detectability durations. Multiply¬ 
ing these times by the respective capture rates, one gets rough 
estimates for the expected average number of LISA sources at 
the Galactic center. However, the MC data permits more precise 
determinations through the following method. Por each event i, 
we computed the evolution of the 20 first harmonics of the rms 
amplitude and selected at each time the one with the highest 
LISA signal-to-noise ratio (SNR). Then we compute 
the time during which there is a harmonic with SNR larger than 
s. Pinally, we get the expected number of sources with a har¬ 
monic component stronger than SNR s by summation, 

^ 1=1 

The summation is realized over all the Nevents events that oc¬ 
curred during some given time interval ATsimui- 

Pig. 3 is the result of this computation. To have reason¬ 
able statistics, we set ATsimui = 10® yrs for MSSs and WDs and 
2 X 10*® yrs for NSs and SBHs. The most striking results con¬ 
cern MSSs. Even though complete tidal disruption has been 
assumed when the particle reaches Rj, the predicted number of 
sources with SNR > 10 is of order 3-5. If sources with SNR 
of 3 can be detected, then LISA should be able to observe of 
order 10 MSSs orbiting the Sgr A* MBH. Note that for such 
long-lived sources, the SNR is proportional to the square root 
of the mission duration, likely to exceed one year. 

4. DISCUSSION 

A major concern regarding this result is the role of tidal in¬ 
teractions which may enter the game before they are strong 
enough to disrupt a star in a single passage. Tides are raised 
at each pericentre passage, converting some orbital energy into 
stellar pulsation. These oscillations dampen into heat through 
mechanisms and on time scales that are still uncertain (McMil¬ 
lan et al. 1987; Kumar & Goodman 1996). Conservatively, we 
assume complete dissipation, thus disregarding the possibility 
of energy being transfered back from the oscillations to the or¬ 
bital motion at next periastron passage, and ineffective radia¬ 
tion of this energy. The ratio of the amount of energy to be 
dissipated (either by gravitational waves or tides) to reach a rel¬ 
ativistic detectable orbit, to the star’s self-binding energy is of 
order (Rs/^p)(c/Vh .)^(1 -e) » 1 where V^, « 600kms“* is the 
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escape velocity from the stellar surface. Consequently, with our 
assumptions, tidal dissipation cannot influence the orbit with¬ 
out first strongly affecting the stellar structure. Hence, to get 
a (over)estimate of the impact of tidal interactions in reducing 
the expected number of detectable stars, we computed for all 
capture events the tidal energy pumped into the star at each pe- 
riastron passage, assuming e ~ 1 and using the n = 1.5 polytrope 
model from McMillan et al. (1987). When the accumulted en¬ 
ergy reached 20 % of the self-binding energy, the star was con¬ 
sidered disrupted. Recomputing with this pessimistic 

procedure, we still get of order 0.5-2 and 4-8 MSS sources 
with SNR >10 and SNR > 3, respectively, see Fig. 3. We con¬ 
clude that tidal interactions can potentialy reduce the number of 
sources significantly but probably not supress them completely. 

We have tested the robustness of our result to several other 
aspects of the modeling, through a series of lower resolution 
simulations (A^p = 10®). Interestingly, our conclusions are not 
affected if we assume a constant stellar formation rate during 
the 10^® years preceding the beginning of the simulation. Ob¬ 
servations of the Arches cluster, close to the Galactic center, 
are hinting to an IMF exponent of a ~ 1.65 for M* > IOMq 
(F iger et al. 1999). If we adopt a= 1.65 for stars more massive 
stars than %Mq, the capture rate is slightly lower. With the ex¬ 
treme assumption that a = 1.65 above Q.SMq, a very SBH-rich 
population is formed. After a few billion years, they dominate 
the central pc from which they have expelled most of the MSSs 
whose capture rate is reduced by a factor of 3. Modifications of 
the low-mass part of the IMF have direct consequences because 
only stars less massive than ~ Q.6Mq may be detected before 
tidal disruption. The IMFs with the largest or lowest proportion 
of low mass stars that are still compatible with local observa¬ 


tions (Kroupa 2001) lead to a number of predicted MSS sources 
higher by 40-80 %, respectively lower by of order 25 %. 

Only low mass MSSs contribute as more massive but less 
denser ones suffer from early tidal disruption. Hence, the MSS 
signal would be too weak to be detected from distant galaxies. 
A MS source would only be marginally detectable during of or¬ 
der 10^ years if located in the Virgo cluster which is too short to 
ensure a significant detection probability. Local group galaxies 
are probably ruled out as well. The BH at the center of M 31 
has a mass of a few Mq so that the frequency of gravita¬ 
tional waves would be too low. Detailed simulations for M 32 
yield a probability for a MS source with SNR above 1 of order 
0.3-0.5 only. If massive BHs are present in dwarf galaxies or 
globular clusters (Gebhardt et al. 2002; Gerssen et al. 2002) and 
their masses follow the well known correlation with the veloc¬ 
ity dispersion of the host galaxy (Tremaine et al. 2002), they 
impose too intense tidal forces in their vicinity. To summarize, 
captured MSSs should be detected at the Galactic center, many 
10 ® years before plunge, but all other sources are predicted to 
be compact remnants in galaxies at distances of a few hundreds 
of Mpc, during the last few months or years of inspiral. 

The draft of this paper has been written during a visit to 
Princeton University Observatory thanks to the hospitality of 
Laurent Eyer. This work was realized under a grant from the 
Swiss National Science Foundation. Complementary support 
from NASA under grant NAG5-10707 is acknowledged. The 
author expresses his gratitude to Kostas Glampedakis who pro¬ 
vided the routines to compute the orbital evolution, to Sterl 
Phinney for interesting discussions and to Scott A. Hughes and 
Tal Alexander for useful comments on the draft of this paper. 
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Fig. 1.— Captures through emission of gravitational waves for our Monte Carlo simulation of the Galactic center, (a) Evolution of the capture rates for the various 
stellar species. Note that only a small number of events have occurred for SBHs or NSs, hence the noisy curves, (b) Orbital parameters at capture for each event 
(which has a statistical weight of 65.5 stars), e is the eccentricity and Rp the pericenter distance (in units of the Schwarzschild radius, Rs = 2.5 x 10“^ pc). Only 
events that occurred in the first lOGyrs of the simulation are included. The surface of disks for MSSs is proportional to the mass of the star. Capture of compact 
remnants are represented with diamonds. 
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Fig. 2.— Gravitational signal for two events from our Sgr A*simulation, a WD and a low-mass MSS. We plot the amplitude versus frequency for the 5 first 

Fourier components of the quadrupolar radiation (Pierro et al. 2001). We use hf where Fobs is the effective duration of observation at frequency /, i.e., 

the mission duration (one year) if the signal is stable over that period. The sharp drop in the signal amplitude at the end of the evolution is due to either the time left 
before swallowing or the characteristic time for frequency increase, /(d//d/)“^, being shorter than the mission duration. The crosses represent the position 1000 
years before plunge through the horizon. Other ticks are for positions 10^, 10^, 10"^, 100, 10, 1 year, 1 month and 1 day before plunge. The segments in dotted line 
for the MSS correspond to a regime where tidal disruption should have occun'ed. The solid V-shaped curve represents LlSA's intrinsic noise (SNR = 1, 

Larson et al. 2000). The dashed line is an estimate of the confusion noise due to unresolved WD binaries in our Galaxy (Bender & Mils 1997). 



Fig. 3.— Expected number of sources of gravitational waves at the Galactic center. This figure shows the number of MSSs, WDs, NSs and SBHs predicted to 
produce a signal above a given signal-to-noise ratio (SNR). For MSSs and WDs, the capture events that occun'ed in the first 10^ years of our Monte Cai'lo simulation 
have been included in the computation (95 and 11 events respectively). For NSs and BHs, we used the first 2 x 10^® years (13 and 22 events). The orbital evolution 
of each captured star, as driven by emission of gravitational radiation around a non-spinning black hole, has been integrated down to plunge instability or tidal 
disruption (Glampedakis et al. 2002) and, at each time, we select the Fourier component of the quadrupolar radiation (Pierro et al. 2001) yielding the highest SNR. 
The upper curve for MSSs is obtained when tidal heating is neglected. The lower curve corresponds to a pessimistic estimate of the decrease in the number of 
sources due to tidal heating, as described in the main text (the number of contributing events is reduced to 27). 






































